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Abstract 

Identified particle spectra of pions, kaons and (anti)protons, and elliptic flow 
Q \ and azimuthal dependence of Bose- Einstein or HBT correlations of identified 

^ ' pions in ^/SNN = 200 GeV Au+Au collisions is analyzed simultaneously using 

^ an ellipsoidally symmetric generalization of the Buda-Lund hydrodynamical 

model. The transverse flow is found to be faster in the reaction plane than 

T 1 

^ i out of plane, which results in a reaction zone that gets slightly more elongated 

' in-plane than out of plane. 

: 1 Introduction 

T— I ■ 

Important information about the properties of extremely hot strongly inter- 
acting matter comes from the observation of azimuthal anisotropies in non- 
^ I central ultra-relativistic nuclear collisions. The second order Fourier com- 

ponent of azimuthal hadron distributions is connected with the azimuthal 
dependence of transverse collective expansion velocity of the bulk matter 
[H |2]. That is in turn determined by the differences of the initial pressure 
gradients in the two perpendicular transverse directions, as well as by the 
initial geometry, the initial velocity and temperature distributions of the fire- 
ball, and the equation of state [21 S] . The anisotropic shape of the fireball 
measured with the help of correlation femtoscopy |S] at the instant of final 
decoupling of hadrons bears information about the total lifespan of the hot 
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matter: with time the originally out-of-reaction-plane shape becomes more 
and more round and may even become in-plane extended [B] . Unfortunately, 
in determining the elliptic flow and azimuthally sensitive correlation radii 
individually two effects — spatial and flow anisotropy — are entangled. For ex- 
ample, the same elliptic flow can be generated with varying flow anisotropy 
strength if the spatial anisotropy is adjusted appropriately [7]. 

In general, the precise way of the interplay between the two anisotropics 
is model dependent. It has been studied and shown to be different within 
the Buda-Lund model [11] than in the Blast Wave model [7]. 

Here we report the results of the Buda-Lund hydro model analyzis of 
azimuthally sensitive Hanbury Brown - Twiss (HBT) radii, using data from 
non-central heavy ion collisions at RHIC Note that the original, axially 
symmetric version of the Buda-Lund model described successfully data from 
central Au+Au collisions at RHIC, as measured by BRAHMS, PHENIX, 
PHOBOS, and STAR collaborations, including identified particle spectra and 
transverse mass dependent HBT radii as well as the pseudorapidity distribu- 
tions of charged particles as first presented in refs. [SI [H] . 

The Buda-Lund model formalism for non-central collisions, including el- 
liptic flow and azimuthal angle dependence of HBT radii has been proposed 
first in [11]. The model is defined with the help of its emission function. In 
order to take into account the effects of resonance decays it uses the core-halo 
model [12]. This ellipsoidal extension of the Buda-Lund model was shown 
before to describe well the transverse mass and the pseudorapidity depen- 
dence of elliptic flow parameter V2 of identified particles at various energies 
and centralities in ref. |10] . 

In the present study, we improve on earlier versions of the Buda-Lund 
model, by scrutinizing the various components using azimuthally sensitive 
HBT data. Eventually we utilize a model that includes as a special case of 
T.S. Biro's axially symmetric and accelerationless exact solution of relativis- 
tic hydrodynamics [13], in contrast to the original, earlier variant, ref. |11] . 
which was based on an ellipsoidally symmetric, but also non-accelerating ex- 
act solution of relativistic hydrodynamics, given by ref. [H]. Similarly to 
ref. [ISj, we utilize here an improved calculation, using the binary source 
formalism, to obtain the observables by using two saddle-points instead of 
only one. This results in an oscillating pre-factor in front of the Gaussian in 
the two-pion correlation function that we take into account for the formulae 
of the HBT radii. Details of the model and the evaluation of the observables 
from it are presented in ref. [22] • 

Azimuthally sensitive HBT radii were also considered recently in cascade 
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models, e.g. in the fast Monte-Carlo model of ref. [T7], or, in the Hadronic 
Resonance Cascade [IB]. 

Data analysis of correlation HBT radii performed earlier with the Blast 
Wave model indicates that the fireball at the freeze-out is elongated slightly 
out of the reaction plane [19], i.e. spatial deformation is similar as in the 
initial state given by the overlap function. This is also supported by the 
theoretical results from hydrodynamic simulations [H1[2U] and URQMD [2T| . 
It sets limitations on the total lifespan. From all previous analyses it seems, 
however, that the final state anisotropy has an interesting non-monotonous 
dependence on collision energy with a minimum at the SPS energies [21]. In 
our analysis of the same data with a different model we observe for the first 
time at RHIC an in-plane elongation of the fireball at freeze-out. 

This presentation is a conference contribution which is based on a more a 
detailed manuscript [22], where we described the azimuthally sensitive Buda- 
Lund model fully and detailed the analytic formulas that were obtained from 
it and were fitted simultaneously to the single particle spectra of identified 
particles, the elliptic flow and the azimuthal angle dependent HBT radii of 
identified pions. 

2 asBuda-Lund hydro compared to data 

Observables like spectra, elliptic flow or Bose-Einstein correlation functions 
have been calculated analytically from the azimuthally sensitive Buda-Lund 
hydrodynamic model, (or asBuda-Lund hydro in a shortened form), us- 
ing a double saddle-point approximation in the integration, as detailed in 
refs. [221 |26]. We have determined [22] the best values of the model pa- 
rameters by fitting these analytic, parametrically given expressions for the 
observables to experimental data with the help of the CERN Minuit fitting 
package. This is possible given that the Buda-Lund hydro model relies on 
exact but parametric solutions of (relativistic) hydrodynamics in certain lim- 
iting cases when these solutions are known, and interpolates between these 
solutions in other cases. So dynamics is mapped to time evolution of model 
parameters, and hadronic observables are sensitive only to the values of these 
model parameters around the time of freeze-out, as was explicitely demon- 
strated recently in ref. [25]. However, penetrating probes, for example the 
direct photon spectra, are known to carry explicit information about the 
equation of state through the cooling history of the continuous radiation of 
these penetrating probes for example direct photons. This property of the 
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model was used recently to determine the equation of state parameters of 
the strongly interacting quark-gluon plasma in 200 GeV Au+Au collisions in 
refs. 1211 Eg 

Data from 20-30% centrality class of 200 AGeV Au+Au collisions pro- 
vided by PHENIX [32l |33] and STAR [SI [35] were used in this analysis. 
The fits were performed simultaneously to azimuthally integrated transverse 
mass spectra of positive and negative pions, kaons, and (anti)protons |32] . 
the transverse momentum dependence of the elliptic flow parameter V2 of 
pions |3l] and to the HBT radii due to pion correlations as functions of 
transverse mass and the azimuthal angle [35]. The results are plotted in 
Figs. 1-3. 

The interpretation of the model parameters is summarized in Table [H 
The two radii Rsx and Rsy correspond to "thermal surface" sizes, correspond- 
ing to distances where the temperature drops to half of its central value, to 
To/2, while parameter Tg corresponds to the temperature of the center after 
most of the particle emission is over (cooling due to evaporation and expan- 
sion). Sudden emission corresponds to the Te = Tq, and Ar — )■ limit. Also 
note that we use /i^, baryochemical potential, calculated from the chemical 
potential of protons and antiprotons as /i^ = 1/2 (/io,p — A^o.p), see Table [H 
The flow profile is linear in both transverse directions but the Hubble con- 
stant is direction dependent, denoted by and Hy in the reaction plane, 
and out of the reaction plane, respectively. 

In Table [21 we present the model parameters obtained from simultaneous 
fits to the data sets. For comparison, results are shown from our earlier 
analysis of 0-30% centrality collisions [36], too, that was performed with a 
previous version of the model corresponding to the axially symmetric limit 
of the current ellipsoidally generalized Buda-Lund hydrodynamic model. 

The general observation is that the Buda-Lund model parameters de- 
scribing the source of non-central reactions are usually slightly smaller than 
those of more central collisions. However, the changes are usually within 2 
standard deviations, therefore the above statement is based on the tendency 
of the parameters, and on some lower energy results not shown here but 
presented in ref. [SB], too. For example, the central temperature in these 
particular non-central reactions is below that of the more central ones. Also, 
the transverse geometrical radii at the mean emission time are considerably 
smaller compared to the more central values. Moreover, the geometric shape 
evolution due to the asymmetric particle transverse flow in-plane (x) an out- 
of-plane (y) directions results in a source more elongated in-plane. Due to the 
smaller longitudinal source size, the parameter corresponding to the forma- 
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tion of hydro dynamic phase is about 10% smaller than that in more central 
collisions, ro(20 — 30%) = 5.4 ± 0.1 fm/c. The elongation in longitudinal 
direction is similarly smaller, A?7(20 — 30%) = 2.5 ± 0.3. In both cases the 
baryochemical potential is found to be small as compared to the proton mass. 
We emphasize again that the observations are based on all the fit results in 
ref. [3S]. 

Note that some of the azimuthally sensitive data have large systematic 
errors that affect the success of fits which we had to take into account. The 
reason for that is the difficulty of precise determination of the event reaction 
plane the data are relative to. Several methods are used by the experiments 
to overcome it and we mention those applied for the selected data. 

The data set we used for fitting V2 was calculated by the four-particle 
cumulants reaction plane determination method that is based on calculations 
of A^-particle correlations and non-fiow effects subtracted to first order when 

is greater than 2. The higher is the more precise the event plane 
determination is, as expected. STAR published two-particle cumulants V2 
data in the same reference, too, but because of the visible deviations between 
the two kinds of data sets and with respect to the comments above we used 
^2(4} data, only. For further details, see ref. [Mj . 

In case of azimuthally sensitive correlation radii, STAR has cast about 
10% possible systematic errors on the data on average. The most likely devi- 
ations were assumed to take effect on the 'side' and 'out' radii of transverse 
momentum of 0.2 GeV/c. The x^/NDF for the full fit, including HBT radii 
with their statistical errors is 269.6/152, which corresponds to a very low 
confidence level. But, when we tested our fits with the above mentioned two 
radii of 'side' and 'out' of transverse momentum of 0.2 GeV/c shifting them 
within their systematic errors (about ± 5%) we could achieve an acceptable 
1% confidence level for the full simultaneous fit. Without the contribution 
of the HBT radii to x^/NDF the confidence level is of an acceptable level of 
5.1%. 

Our results were compared also to results of refs. [281 E] of the azimuthally 
sensitive extension of the Blast Wave model in the detailed write-up, ref. [22] . 

3 Conclusions and outlook 

The ellipsoidally symmetric generalization of the Buda-Lund hydrodynamic 
model compares favourably to the identified particle spectra and to the el- 
liptic fiow and azimuthally sensitive Bose- Einstein correlation radii of iden- 
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Figure 1: Buda-Lund model fits to ^/snn = 200 GeV Au+Au data of ref. [32] 
on azimutlially integrated transverse momentum spectra of negatively (top 
figure) and positively (bottom plot) charged particles. 
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Figure 2: Buda-Lund model fit to RHIC 200GeV Au+Au data on f 2 elliptic 
flow of pions, data from ref. |34j . 



Buda-Lund 


parameter description 


To 


Temperature in the center at tq 


Te 


Temperature in the center at Tq + At 




Baryochemical pontential in the center at tq 


Rx 


Geometrical size in direction x 


Ry 


Geometrical size in direction y 


Rxs 


Thermal size, where T = To/2 in direction x 


Rys 


Thermal size, where T = To/2 in direction y 


Hx 


Hubble constant in direction x 


Hy 


Hubble constant in direction y 




Mean freeze-out proper-time 


At 


Distribution width in proper-time r 


At] 


Distribution width in space-time rapidity 77 




Baryochemical potential. 



Table 1: Description of the parameters of the asBuda-Lund hydro model. 
The baryochemical potential is evaluated as = 1/2 (/^o,p — /^o,p)- For 
details see ref. 1221. 
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Figure 3: Buda-Lund model fits to RHIC 200GeV Au+Au data of ref. [35] 
on HBT radii as functions of transverse mass, for different azimuthal angles 
(top). The same data and the same fit is also shown as a function of the 
azimuthal angle for different values of the transverse mass (bottom). 
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tified pions. From model fits to 20-30 % central Au+Au collision data at 
y^s/vjv = 200 GeV at mid-rapidity, the source parameters characterizing these 
non-central ultra-relativistic heavy ion reactions were extracted. 

The results of our analysis indicate that the central temperature Tq = 
174 ± 6 (stat) MeV in the 20-30% centrality class is somewhat lower than 
that in more central collisions, where an earlier analysis found Tq = 196 ± 
13 (stat) MeV. We have found that the transverse flow is stronger in the 
reaction plane than out of plane with Hubble constants = 0.158 ± 0.002 
and Hy = 0.118 ± 0.002. The almond shape of the reaction zone initially 
elongated out of plane gets slightly elongated in the direction of the impact 
parameter by the time the particle emission rate reaches its maximum. The 
effect is reflected by the geometrical radii in the two perpendicular directions 
at that time, Rx{iTi — plane) = 9.5 ± 0.5 (stat) fm, i?j^(out — plane) = 7.0 ± 
0.2 (stat) fm. As far as we know, this study is the first one where an in-plane 
extended source has been reconstructed from simultaneous and reasonably 
successful hydro model fits to identified particle spectra, elliptic flow and 
azimuthally sensitive HBT data in 200 AGeV Au+Au collisions at RHIC. 
It is a remarkable property of this hadronic final state, that the ratio of the 
Hubble constants is approximately the same as the ratio of the geometric 
source sizes: H^/Hy = Rx/Ry within the errors of the analysis. 
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Buda-Lund 
parameters 


Au+Au 200 GeV 
central (0-30%) 


Au+Au 200 GeV 
non-central (20-30%) 


To [MeV] 
Te [MeV] 
f^B [MeV] 


196 ± 13 
117 ± 11 
31 ± 28 


174 ± 6 
130 ± 6 
27 ± 16 


Rx [fm] 
Ry [fm] 
Rsx [fm] 
Rsy [fm] 


13.5 ± 1.7 

= Rx 

12.4 ± 1.6 

Rsx 


9.5 ± 0.5 
7.0 ± 0.2 

12.8 ± 0.8 

16.9 ± 1.6 


Hx 

Hy 


0.119 ± 0.020 

Hx 


0.158 ± 0.002 
0.118 ± 0.002 


To [fm/c] 
At [fm/c] 
Ar/ 


5.8 ± 0.3 
0.9 ± 1.2 
3.1 ± 0.1 


5.4 ± 0.1 

2.5 ± 0.2 
2.5 ± 0.3 


XVNDF 


114/208 = 0.55 


269.4/152 = 1.77 



Table 2: Source parameters from simultaneous fits to PHENIX and STAR 
data of Au+Au collisions at a/J/vtv = 200 GeV, as given in Figs. 1-3, obtained 
by the Buda-Lund model. For non-central data the value of x^/NDF refers 
to fits with statistical errors, only. 
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